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I discuss a QCD sum rule determination of the four-quark operator matrix elements contributing to the Aj, 
inclusive decay rates at C(m^~ 3 ). The results suggest that l/m 3 corrections are not responsible of the observed 
difference between the lifetimes of Aj, and B^. 



1 Introduction 

The difference between the measured lifetimes of 
the A;, baryon and Bd meson EJ: 



r(Ag)/r(B°) = 0.78 ± 0.04 



(1) 



represents an intriguing problem of the present- 
day heavy quark physics. Indeed, the 20% differ- 
ence contradicts the expectation that, at the scale 
of the b quark mass, the spectator model should 
describe rather accurately the decays of hadrons 
Hq containing one heavy quark. 

A calculation of the ratio t(Ai ) )/t(B c [) can be 
attempted using a field theoretical approach de- 
veloped for the analysis of the inclusive weak de- 
cays of the hadrons Hq a. The method is based 
on the expansion in powers of rag 1 , in the frame- 
work of the Wilson OPE. The widths are expressed 
in terms of hadronic matrix elements of high di- 
mensional quark and gluon operators; since such 
matrix elements can be responsible of the large 
difference between t(A&) and r(Bd), it is interest- 
ing to compute them by nonperturbative methods 
such as QCD sum rules. 

The main aspects of the QCD analysis of the 
inclusive decays of hadrons Hq can be summarized 
considering the transition operator T(Q Xf — * 



f = i d^x T[C w (x)C\ v (0)] 



(2) 



which describes an amplitude with the heavy 
quark Q having the same momentum in the ini- 
tial and final state. Cw is the effective weak La- 
grangian governing the decay Q —y Xf. The in- 
clusive width Hq —y Xf can be obtained from 



X, 



21m <Hq\T\H q > 



(3) 



The large energy release in the heavy quark de- 
cay permits an expansion of T in terms of local 
operators Oc 



T 



(4) 



with Oi ordered according to the dimension, and 
the coefficients Cj containing appropriate inverse 
powers of the heavy quark mass tuq. The lowest 
dimension operator in ([|) is O3 = QQ. The next 
gauge and Lorentz invariant operator is the D = 
5 chromomagnetic operator Oq — Q^cr^G^Q, 
whose hadronic matrix element 



V 2 g(Hq) 



< H Q \Q^G^Q\H Q > 



2M 



(5) 



measures the coupling of the heavy quark spin to 
the spin of the light degrees of freedom in Hq , and 
therefore is responsible of the mass splitting be- 
tween hadrons belonging to the same se multiplet 
(s£ is the total angular momentum of the light de- 
grees of freedom in Hq ) . In the case of b- flavoured 
hadrons this mass difference has been measured, 
both for mesonsL(Af B . - M B = 45.7 ± 0.4 MeV) 
and S 6 baryonsi (M E; - M s „ = 56 ± 16 MeV). 

The matrix element of QQ over Hq can be 
obtained using the following expansion stemming 
from the heavy quark equation of motion: 



QQ = Qj°Q 



O 



G 



2ml 



2m 2 Q 



0(m Q 3 ) 



(6) 



O n is the kinetic energy operator T: 
whose matrix element 



Q(iD) 2 Q 



2 <Hq\Q( 1 DYQ\H q > 

VAHq) = rTTT (') 

2Ad H Q 

measures the average squared momentum of the 
heavy quark inside Hq. 
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At the order C(toq ) four-quark operators ap- 
pear in (Q): 

0% = QTq qTQ (8) 

with r an appropriate combination of Dirac and 
color matrices. 

The resulting expression for the width 
T(H Q -> X f ) reads: 



r(H Q ^x f ) = r f \aZ + ^. + \. 



(9) 



A{ and Tg depend on the final state Xf, A{ in- 
clude perturbative short-distance coefficients and 
nonpcrturbative hadronic matrix elements incor- 
porating the long range dynamics. The lead- 
ing term in (§) corresponds to the partonic pre- 
diction TP art (H Q X f ) = T f Q A f Q , with A f = 
1 + cJ a s /-K + 0{a 2 s ) and Tq oc uiq; differences 
among the widths of the hadrons Hq emerge at 
the next to leading order in 1/mg, and are re- 
lated to the different value of the matrix elements 
of the operators Oi of dimension larger than three. 

It is important to notice the absence of the 
first order term tUq 1 in (^H 

The D — 5 operators Oq and O n are SU(3) 
singlets; on the contrary, the D — 6 operators in 
(H) give different contributions when averaged over 
hadrons belonging to the same SU (3) light flavour 
multiplet, and therefore they are responsible of the 
different lifetimes of, e.g., B~ and B s , A b and 3(,. 
The spectator flavour dependence is related to the 
mechanisms of weak scattering and Pauli interfer- 
ence □, both suppressed by the factor tUq 3 with 
respect to the parton decay rate. 

As for differences between mesons and 
baryons, they could already arise at C(mg 2 ), due 
both to the chromomagnetic contribution and to 
the kinetic energy term in (||). In particular, the 
kinetic energy term is responsible of the difference 
for systems where the chromomagnetic contribu- 
tion vanishes, namely in the case of A& and Sf, 
having the light degrees of freedom in S— wave. 
However, the results of a calculation of /z^ fox 
mesons □ and baryons □ support the conjecture 
that the kinetic energy operator has the same ma- 
trix element when computed on such hadronic sys- 
tems. □ The approximate equality of the kinetic 
energy operator on Bd and Af, can also be inferred 
by using mass relations fJ^(Ab) — ^(Bd) — 
0.002 ± 0.024 GeV 2 . Then, differences between 



meson and baryon lifetimes should occur at the 
rag 3 level, thus involving the four-quark operators 
in (||). They can be classified as follows l!3: 



U V-A 
U S~P 
rpq 

1 V-A 



rpq 

± S-P 



= QnqL q~LQn 

= Qhi^qL qhi^Qh 

n Xa - X "n 
= Qii^qL qh-^Qn 



(10) 



(qn,i 



_ 1±75 



q, X a Gell-Mann matrices). 



For mesons, the matrix elements of the oper- 
ators in ([To]) can be computed by vacuum satura- 
tion approximation: 



< B q \O q v _ A \B v > VS A 

< B q \T v _ A \B q >vsa 







- (11) 
(12) 



etc. Such an approximation cannot be employed 
for baryons, where a direct calculation is required. 

For A;,, a simplification can be obtained IMj in- 
troducing 



Tl qli'qi 



(13) 



and 



U S-P 



QW R qlQ l R (14) 

(i and j color indices): the Af, matrix elements of 
the operators in (10) can be expressed in terms 



of < A b \Olr_ A \A b > and < A b |0^_ A |Af, >, mod- 
ulo 1/toq corrections contributing to subleading 
terms in the expression for the inclusive widths. 
Parametrizing the matrix elements 



(O q V_ A )A b 



< A b \O q v _ A \A b > 



PbMb, 



2M Ab 48 
< A 6 |0«_JA 6 > = -B < A b \6 q v _ A \A b > 



(15) 



corresponds to {Oy_ A ) A 



(16) 

200 MeV and r = 
q \- - 4.4 x 



one has that, for fg 
1, eq-dH 

10" 3 GeV 3 . Quark models predict {6 q v _ A ) Ac ^ 
0.75 — 2.ELx 10 -3 GeV 3 , corresponding to_r ~ 
0.2 — 0.6EJ; larger values can be obtained E3 us- 
ing the mass splitting E£ — Ef, and E* — E c : 
r ~ 0.9 — 1.8. In the valence quark approxima- 
tion B = 1. 

A value of r: r ~ 4 — 5 woukLexplain the dif- 
ference between r(Af,) and r(Bd)i3. A calculation 
by QCD sum rules, however, seems to exclude this 
possibility. 
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2 {O v _ A )\ b from QCD sum rules 

An estimate of (Oy_ A ) Ab can be obtained by QCD 
sum rules, in HQET, analyzing the three-point 
correlator 



W cd {lo,uj') = (1 + ^) cd \\{lo,lo') 



dxdy <Q\T[J c {x)O q v _ A , 



iujv-x—iuj v-y 



(17) 



of the baryonic currents J{x) and J(y) and of 
the operator O v _ A in (p^|). The variables w, a;' 
are related to the residual momentum of the cur- 
rents: — mijV^ + lov^ . Choosing J with a non- 
vanishing projection on the A& state 



<0\J c \At(v) >=/a 6 (^)c 



(18) 



(ip v is a spinor for a A& of four- velocity v), the pa- 
rameter /a 6 representing the coupling of J to A&, 
the matrix element (O v _ A )\ b can be obtained by 
saturating the correlator (^) with baryonic states, 
and considering the double pole contribution in 
the variables u> and u>': 



Il h (uj,uj') 



Pv-a)a, 



2 (A Ab -w)(A/ 



at the value lu = u>' = A/ 



(19) 

The mass parameter 



A\ b represents the binding energy of the light A& 
degrees of freedom in the static color field gener- 
ated by the 6— quark: M\ b = nit, + A\ b ; it must 
also be derived by QCD sum rules. _ . 

A possible interpolating field for A&llia reads: 

Jc(x) = e» k (q Ti (x)Trq^x))(h k v )c(x) (20) 

where T means transpose, i, j and k are color in- 
dices, and C is the Dirac index of the effective 
heavy quark field h v (x) ; r is a light flavour matrix 
corresponding to zero isospin. In the — > oo 
limit, the light diquark in the Af, is in a relative 
+ spin-parity state; therefore one can use in (Uw: 



C 75 (l + M ) 



(21) 



C being the charge conjugation matrix and b a 
parameter. Arguments .can be found in favour of 
the choice b = 1 in (§|)B. pn 
The coupling /a 6 can be computed fromHH-3 

H CD {u) = {l+f) CD H{Lu) 
= i [ dx< Q\T[Jo(x)J D (p)]\0 > e iuv - x (22) 



In the Euclidean region the correlation functions 
( |i"7j ) and (^2|) can be computed in QCD, in terms 
of a perturbative contribution and vacuum con- 
densates: 

*° PB M = [ dada' ""frfl (23) 
J (a - u)){a' - u/J 



H ope (uj) = / dap 
J [c 



■to) 

with the spectral functions given by : 



(24) 



pn(a,a>) = p% ert \a, a') + p^{a, a') < qq > 

(£>=4)/ ,* ^ f* s 2 

+ Pn i^ cr )< —G > 

+ pg= 5 \a,a')<qg*Gq> 

+ p ( °= 6 \a,a')(<qq>f + ... (25) 

and a a similar expression for pnio). 

The various terms in pu(a,cr') and ph(ct) can 
be found intllia. 

A sum rule for (O v _ A )A b can be derived by 
equating the hadronic and the OPE representa- 
tions of the correlator (|lj), and modeling the con- 
tribution of the higher resonances and of the con- 
tinuum in IP as the QCD term outside the region 
< u> < uj Ci < oj' < u) c , lu c being an effec- 
tive threshold. After a double Borel transform one 
gets: 



lk {l + b f e -^ { o v _ A ) hb (26) 



dcrda'e ^ p n (cr,cr') (27) 



where E is a Borel parameter. The threshold co c 
can be fixed in the QCD sum rule determination of 
f\ b and Aa 6 □: lo c = 1.1 — 1.3 GeV; in correspon- 
dence one obtains f Ab = (2.9 ± 0.5) x 10~ 2 GeV 3 
and A Ab = 0.9 ±0.1 GeV. 

The resulting sum rule is depicted in fig.l. A 
stability window is observed for E > 0.2 GeV. In 
the duality region E ~ 0.2 — 0.3 GeV (the same 
region considered in the QCD sum rule analysis of 
/a„ and /4(A h )Q) we find: 



(Oy_ A ) Ab ~ (0.4 - 1.20) x 10~ 3 GeV 3 



(28) 



which corresponds to r ~ 0.1 — 0.3. This result is 
confirmed by an analysis based on the assumption 
of local quark- hadron duality, that amounts to cal- 



culate the matrix elements of (O 



v-A/M 



and / A 
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0.1 0.150.20.250.30.350.40.450.5 



E (GeV) 



Figure 1: Sum rule for the matrix element (6ui}a 4 



by free quark states produced and annihilated by 
the baryonic currents in (17) and (|2l|L_and then 
averaging on a duality interval in w,u/li3. 

Considering, finally, the parameter B in 
eq.([l6|), one gets B = 1 since, in this computa- 
tional scheme, only valence quark processes are 
taken into account. 



3 Conclusions 

Within the uncertainties of the method, QCD sum 
rules predict small values for the matrix elements 



V-AI 



and 



i)A(,) comparable with the 



outcome of constituent quark models. The conclu- 
sion is that the inclusion of 1 /nig terms in the ex- 
pression of the inclusive widths does not solve the 
puzzle represented by the difference between r(Af,) 
and T(Bd): using the formulae inlij for the lifetime 
ratio, the value in (Ea) together with B = 1 gives: 



(29) 



-(A b )/ T (B d ) > 0.94 



It seems unlikely that order tuq terms can 
solve the problem. We must conclude that, if 
the measurement of r(Ab) and r(Bd) will be con- 
firmed, a reanalysis of the problem will be re- 
quired. In particular, it has been proposed to 
consider the possibility of the failure of the as- 
sumption, made in the calculation of the lifetimes, 
of local quark-_haihpn duality in nonleptonic in- 
clusive decays Ej'EZnlfl. Meanwhile, it is interest- 
ing that new data are now ^available for other 
b— flavoured hadrons, e.g. S^El, although with er- 
rors too large to perform a meaningful comparison 
with At,. Such new information will be of great im- 
portance for the full understanding of the problem 
of the beauty hadron lifetimes. 
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